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We discuss two related aspects of the lattice version of the heavy quark effective theory (HQET). They are 
the effects of heavy quark modes with large momenta, near the boundary of the Brillouin zone, and the 
renormalization of the lattice HQET. We argue that even though large momentum modes are present, their 
contributions to heavy-light bound states and perturbative loop integrals are dynamically suppressed and vanish 
in the continuum limit. We also discuss a new feature of the renormalization of the lattice HQET not present 
in the continuum theory, namely that the classical velocity is finitely renormalized. 



1. INTRODUCTION 

The Isgur-Wise form of the heavy quark effective 
theory (HQET) is based on the separation of mass 
scales and on the physical idea that the motion of a 
very heavy object is only slightly perturbed by its 
interaction with much lighter objects. The idea of 
the Isgur-Wise limit is that a heavy quark's 
momentum can be usefully decomposed as 
= Mv^ + , where v is a constant but arbitrary 

four- velocity normalized to = 1 .''^ In the 
M->oo limit, the fermion propagator simplifies to 

1 ^ 1^-^ (1) 
yP-M 2vp 

The power of this limit is that the large momentum 
Mv disappears, leaving behind only terms involving 
V and p, which are o(l) in the limit. Another 
important simpUfication is that the Dirac matrix 
structure of the propagator becomes trivial. This is 
responsible for the spin symmetry that appears in the 
Isgur-Wise limit. While the same physics lies behind 
both the static quark analysis^ widely used in lattice 
simulations and the Isgur-Wise form of the HQET, 
the fact that v can be non-zero in the Isgur-Wise 
form allows the calculation of processes with large 
momentum transfer, on the order of the heavy quark 
mass, such as the Isgur-Wise universal form factor 



describing the semileptonic decay of one particle 
containing a single heavy quark into another. 

We have given earlier a lattice implementation of 
the HQET.* It shares with the continuum version the 
possibility that v can take any value, and yet the 
momenta that enter simulations or perturbative loops 
remain small, as is appropriate for lattice analysis. 

This lattice heavy quark effective theory (LHQET) 
has some unfamiliar features which we will discuss 
in this talk. One is the that there are modes of the 
heavy quark field that grow rather than decay with 
Euclidean time. A related property of the LHQET is 
that the heavy fermion modes are doubled. We 
will argue that the potential unphysical effects of 
these properties are absent, and that the same 
physical mechanism operates in both cases. 

The other feature of the lattice HQET we will 
address here is that the classical velocity v is 
(finitely) renormalized. We will show how this new 
renormalization arises from the lattice discretization 
of the reduced Dirac equation for the heavy quark 
propagator. We will give a one-loop perturbative 
estimate of its magnitude, and also compute it in a 
lattice simulation. The estimates are in fairly good 
quantitative agreement. 

2. THE LATTICE HQET 

The heavy quark effective theory is efficiently 
formulated' by factoring the M^oo singular behavior 



from the field of a heavy quark, leaving a reduced 
operator . 

h'^\x) = e '"^"-^ J_^rv (2) 

There is an independent field for each classical 
velocity. The Lagrangian for h is 

= h'-"\x) ivD h^'Xx) 

and its propagator satisfies the reduced Dirac 
equation 

-ivDS = 5 (4) 

Here D is the covariant derivative. The 
discretization by which we defined the lattice 
HQET is 

Vg [ U(x,x+i) S{x+i,y)- S{x,y) ] 

-U(x,x-(i) S(x-(i,y) ] 
- 8(x,y) 

The use of an asymmetric forward time difference 
facilitates implementing the requirement that heavy 
quarks propagate only forward in time. A 
symmetrical difference could also be used, at the cost 
of a small increase in computation. 

The free inverse propagator corresponding to the 
forward-time symmetrical- space discretization is 

S^^'-'ip) - v,(e + E v^sinp^ (6) 

n-i 

The balance of this talk will be devoted to a study 
of some rather unfamiliar structural aspects of the 
theory defined by this reduced discrete Dirac 
equation. 

3. THE LATTICE HQET SPECTRUM 

We will discuss together fermion doubling and the 
asymptotic behavior of the reduced heavy quark 
propagator. Doubling occurs because the 
symmetrical first spacial difference results in the 



presence of modes with large residual momentum but 
small residual energy. Growing modes occur 
because the reduced heavy quark Euclidean-space 
propagator is obtained from the full propagator by 
factoring out the zero residual momentum decay. 
This results in those modes whose residual momenta 
are opposite to the classical velocity growing with 
Euclidean time. In the heavy quark limit, the rate of 
decay of the mode with residual momentum p is 
given by the finite residual energy 

]lM^+iMv+pf- m/i + v' ^ (7) 

When this is negative, the "falloff becomes growth. 
Conversely, at fixed Euclidean time, the heavy quark 
propagator grows with momentum in directions 
opposite to V. 

The growing modes do not lead to physical effects 
or loop divergences in the LHQET because they are 
always accompanied by modes of either a light 
Wilson quark or a gluon whose momentum is at the 
edge of the Brillouin zone. The suppression of 
divergences associated with the growing heavy quark 
modes is seen explicitly in perturbative loop 
diagrams.* It occurs because the rate of falloff of 

the modes of the light quantum is / 1 p\^ + m^ , while 
according to Eq. (7), the rate of heavy quark 
propagator's growth is only a finite fraction of | p 
(since |v|/Vp is strictly less than 1). In the 

Euclidean energy plane, these effects appear as a 
result of the integration contour, determined by the 
heavy quark's propagating only forward in time. 

The absence of physical manifestations of the 
heavy quark doubler modes is illustrated by lattice 
free field theory.' The falloff of a reduced 
composite propagator for a "meson" made of a free 
heavy quark and a free light quark with momentum 
P - Mv+ p is given by the minimum of the 
combined energies of a light quark with momentum 
li and a heavy quark with momentum Mv + p- Ic, 
reduced by the heavy quark energy for vanishing 
residual momentum, Mv^. In two dimensions the 

minimization can be carried out analytically. The 
minimum is unique; only one heavy quark mode is 
selected. In the heavy quark limit, the minimum 



occurs at 



and is the correct residual 



energy for an object with mass (M + m), classical 
velocity v, and residual momentum p^, namely 



(Mvj+ Pj)- + (M + mf - Mv^ 



(8) 



4. RENORMALIZATION OF THE CLASSICAL 
VELOCITY 

A new feature of the HQET on the lattice is the 
possibility that the classical velocity is renormalized. 
This is a finite lattice artifact directly related to the 
reduction of the spatial symmetry from (9(4) to the 
hypercubic group and can be estimated both in lattice 
perturbation theory and in simulations. The full 
inverse propagator for a heavy quark with classical 
velocity v is given by 



s -\p) - s^'^-\p) - Up) 



(9) 



4-vectors, such as v^^ , which have exactly the same 
lattice transformation properties as v^^. Such terms 
are present in even the simplest 1 loop contribution 
to E, and they cause the classical velocity to be 
shifted by a finite amount. 

We may use the change in energy for small 
residual momenta (Eq. (7)) to define the physical 
(renormalized) classical velocity without resort to 
perturbation theory. Of course, the quark propagator 
is gauge dependent, and so it is not an appropriate 
quantity to directly simulate. However, the 
propagator of a gauge invariant composite heavy- 
light meson propagator has the same dependence on 
the classical velocity as its heavy quark component, 
and is suitable for direct simulation. The falloff of 
a heavy-Ught propagator is controlled by the location 
of the closest zero (in the residual energy plane) of 
the inverse propagator. From Eq. (1) or (6) this 
gives 



V 



iphys) 



dEiv^°\p) 



(12) 



p = 



where in perturbation theory E is the sum of all 
proper self mass diagrams. Denoting the derivative 
of E with respect to p at zero by X^, the 
renormalized classical velocity is given by 



= V - X 



(10) 



"7 -1 .,, ren 



The requirement that the physical classical velocity 
be normalized to 1 fixes both the renormalized 
classical velocity and the wave function 
renormalization constant Z. 



Z = 



1 



(V, 



(0) 



X y 



(11) 



Z (v'"> - X ) 



We have estimated the magnitude of the classical 
velocity renormalization both perturbatively and by 
simulating Eq. (12). The perturbative calculation 
follows the method of Ref. 6, with the replacement 
of the backwards time derivative used there with a 
forward derivative. For the simulation, we used a 
preUminary subset of 24^ x 48 lattices with P = 6.1 
and K = .154 generously provided by the Fermilab 
Collaboration.** In carrying out the simulation, it 
proved useful to develop the propagator in an 
expansion in powers of v. 



V /v 

V.I v^. 



M{t,p) = v7' v^' ^r^' M{t,p,m) (13) 

This is effectively a transverse hopping expansion, 
because computing the propagator nij units from the 
origin in the i* direction introduces factors 
of V. . For fixed v and p, the asymptotic decay of 

the propagator is 



In the continuum, 0(4) invariance forces X to be 
proportional to v, the only available 4-vector, and so 
its only effect is to shift Z from 1. On the lattice, 
however, there are additional linearly independent 



M{t,p) ~ C(p)e -''^P'" 



(14) 



The normalization and energy are also developed in 
analogous series. Parity invariance implies that 



terms even (odd) in the i* momentum component 
come only from power series coefficients with even 
(odd) values of nij. The logarithmic derivative of the 
asymptotic propagator with respect to the residual 
momentum at the origin identifies, through Eq. (12), 
the physical classical velocity as the coefficient of t. 
The statistical signal seems to deteriorate quickly as 
higher order expansion coefficients are considered, so 
we restricted our preliminary simulation analysis to 
first order in the bare classical velocity. This gives 



dM(t,p,m. = 1)1 dp. 



p = o 



M(f,p =0,w, =0) 
dE{p,m. = 1) 



(15) 



dp, 



t + const 



p=0 



The momentum derivative was approximated by the 
symmetrical first difference, and the ratio of = 1 
to m,- = propagators was fit out to 6 units of 
Euclidean time separation. 

Both the simulation and perturbation theory 
indicate that the physical classical velocity is smaller 
than its input, or bare value. The comparison of the 
shifts in the classical velocity, using the mean-field 
tadpole improved value for the coupling' 
corresponding to P = 6.1 in the perturbative 
evaluation is as follows: 



{ren) 



1 loop 



iphys) 



{input) 



.72 



.6 + .1 



Jsimul 



These values are quite consistent, but perhaps 
surprisingly far from unity. 



arguments, we argued that neither the asymptotic 
growth of some modes of the heavy quark 
propagator, nor the naive appearance of heavy 
fermion doubler modes leads to physical 
consequences. 

We also quantitatively examined the (finite) 
renormalization of the classical velocity. This 
renormalization is peculiar to the LHQET, and does 
not occur in the continuum. At the values of the 
lattice coupling at which QCD simulations are 
typically performed today, this seems, on the basis of 
this very preliminary analysis, to be a significant 
effect. We of course plan to refine this analysis in 
the near future. The renormalization of the classical 
velocity must be taken into account in any attempt to 
use the LHQET to compute the Isgur-Wise function 
or other classical-velocity-dependent quantities, 
strictly in the heavy quark limit, on the lattice. As 
one evident example, the classical velocity rescaling 
changes the inferred values of the slope of the Isgur- 
Wise function at the origin. 
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5. CONCLUSIONS 



We have discussed in this talk some characteristics 

of the lattice implementation of the heavy quark 
effective theory. By means of essentially kinematic 



